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Abstract. We present photometry, photometric redshifts and extra galactic 
number counts for uhra deep 15 micron mapping of the gravitational lensing 
cluster Abell 2218 (A2218), which is the deepest image taken by any facility 
at this wavelength. This data resolves the cosmic infrared background (CIRB) 
beyond the 80% that blank field AKARI surveys aim to achieve. To gain an 
understanding of galaxy formation and evolution over the age of the Universe a 
necessary step is to fully resolve the CIRB, which represents the dust-shrouded 
cosmic star formation history. Observing through A2218 gives magnifications of 
up to a factor of 10, thus allowing the sampling of a more representative spread 
of high redshift galaxies, which comprise the bulk of the CIRB. 19 pointed 
observations were taken by AKARPs IRC MIR-L channel, and a final combined 
image with an area of 122.3 square arcminutes and effective integration time of 
8460 seconds was achieved. The Scr sensitivity limit is estimated at 41.7 /xJy. 
An initial 5cr catalogue of 565 sources was extracted giving 39 beams per source, 
which shows the image is confusion limited. Our 15 micron number counts 
show strong evolution consistent with galaxy evolution models that incorporate 
downsizing in star formation. 



1. Data Reduction 

Our data consists o f 19 pointed observ ation of A2218, t aken with the 15 micro n 
filter of the IRC-L (lOnaka et alJl2nn7l ) aboard AKARI (|Murakami et alJ 120071 ). 



The data was reduced using the standard IRC pipeline, version 20070912. The 
pipeline's median sky subtraction was utilized but resulted in dark areas, signif- 
icantly around the brighter sources. To deal with this issue, and the remnants 
of scattered light persisting post-pipeline, a further median sky subtraction of 
the background areas was performed. Hot pixels were masked and removed, and 
the remaining bad pixels were addressed using an IDL sigma filtering routine. 
Figure 1 shows four corresponding postage stamp sections taken from frame 19 
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illustrating the post-pipeline output, a median filtered mask, the median sub- 
tracted result and the sigma filtered result. The resulting 19 images were aligned 
with Aladin and IDL's HASTROM, and then average combined to give the final 
L15 image. 




Figure 1. Postage stamp image sections showing, from left to right, an 
example of the post pipeline output, median filtered mask, median subtracted 
image and sigma filtered image. 



1.1. Source Extraction and Completeness 

A 5a source extraction was performed with DAOFIND (jStetsonl The 



resulting catalogue was eyeballed to eliminate any spurious detections, giving a 
final number of 565 sources. A Monte Carlo completeness test was performed 
using an IDL routine written to convolve a normalized empirical PSF with the 
final L15 image, to create artificial sources at random positions. The artificial 
sources were placed sufficiently apart from one another and known sources to 
avoid self-confusion. The test was performed for 80 flux bins, covering the range 
of flux densities for detected sources within the L15 image. In order to reduce 
statistical errors the test was repeated until an effective 18452 sources per bin 
was achieved. The completeness test results show the L15 image is 50% complete 
to 30.7 /iJy and 80% complete to 39.4 /iJy. 

1.2. Sensitivity 

Aperture photometry was taken at random positions on the flnal image, exclud- 
ing the edges and the brightest source, and the results were plotted as a his- 
togram of flux density against number, see Figure 2. The resulting distribution 
has an asymmetric tail that signifies the contributions from bright sources. The 
combined confusion and detector noise can be represented by fitting a Gaussian 
to the histogram, with a standard deviation of 8.33 //Jy giving a 5a sensitivity 
estimate of 41.67 //Jy. The mean rms per pixel of the L15 noise map is 2.20 fiJy. 




Figure 2. Distribution of random aperture photometry of the L15 image, 
fitted with a Gaussian of standard deviation 8.33 /iJy. 



Deepest image at a wavelength of 15 microns 
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2. Band-merged catalogue 

Multi waveband data of A2218 was used to identify counterparts of the L15 
5cr source catalogue. HST WFP2 F450, F606 and F814, Palomar 200inch Hale 
U, V, b, i and INGRID WFC Ks and J images of A2218 were provided by 
Ian Smail. Spitzer IRAC Ch 1 to 4 data was obtained via Leopard and com- 
bined with Mopex. An AKARI Sll image was provided by Myungshin Im and 
Jongwan Ko, and a Spitzer MIPS 24 /xm image was provided by Eiichi Egami. 
The counterparting procedure identified an extra 368 sources and 3 spurious 5a" 
detections, giving a combined counterpart catalogue of 930 sources. 



3. Photometry 



Aper ture photometry of the ba source catalogue was taken with PHOT (jStetson 
19871 ^ using an aperture of 5.96" and a sky annulus of radii 19.07" and 31.0". An 
aperture correction of 1.30, derived using a growth c urve correction met hod, was 
applied and the IRC data user's manual version 1.4 ( Lorente et al. 20081 ) conver- 
sion factor of 1.691 was used to convert from ADU to /jJ y. For the HST ima ges, 
photometry was obtained from the published catalogue o f ISmail et al.1 (|200lh via 
the NED database. The remaining images were subject to a growth curve cor- 
rection method to obtain corrected aperture photometry, using a routine written 
in IDL. For each image one or two mean growth curves were empirically con- 
structed. These curves were used to calculate aperture radius and correction for 
each remaining source. IDL's APER was used to take the subsequent aperture 
photometry. Our Ks band photometr y was compared, w here available, to previ- 
ously published Ks band photometrv lSmail et al. ( 200ll ) and showed a less than 
2% difference. 



4. Photometric Redshifts 



We used EaZy (jBrammer. van Dokkum CoPDill2008l ') to gain photometric red- 
shift estimates for our 5o" source catalogue. EaZy utilizes a minimum SED 
fitting method, which is suitable for data sets with no available spectroscopic 
redshifts (Zspec) or a biased set of Zspec. In our case the majority of the small 
Zspec available are biased at the cluster distance. The theoretical SED tem- 
plates, used by EaZy, are based on semi- analytical models, and a linear combi- 
nation of templates can be fitted simultaneously. EaZy gives the option to apply 
priors, aimed at breaking the template colour degeneracies seen with increasing 
redshift. A comparison of the resulting redshift estimates for our sources with 
known Zspec shows that applying priors gives an improved correlation of ap- 
proximately 10%. Our spectra we r e also fitted using the photometric redshift 
code illustrated in iNegrello et al.l ( 20091 ). Photz from here on. This code is 



uniquely optimised for fitting mid-to-far -infrared PAH and silicate features seen 
i n starburst SEDs. Starburst templa te ( Takagi et al.ll200"^ ) and AGN template 



( Efstathiou &: Rowan-Robinsonll 19951 ) components were simultaneously fitted by 



photz. A comparison of the Photz and EaZy redshifts estimates for sources with 
prominent mid-to-far-infrared features shows a correlation of around 0.8. 
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Number Counts 



Differential number counts (dN/dS), corrected for incompleteness, were taken 
for our 5(7 catalogue and normalized to a Euclidean slope. These counts were 
corrected for flu x amplification by applying magnification factors obtained via 
LENSTOOL (Julio et alJ [20071 ^ . Figure 3 compares our corrected and uncor- 
rected counts with previously published differential number counts. Our counts 
corrected for lensing show an upturn arou nd 2 m Jy and peak arou nd 0.4 mJy, in 
agreement with prev i ous c ounts (e.g.. .Elbaz et al.lll999l ) and the iPearson et al" 
(|2007l 1: iPearson. C. I (|2009l ') model. 




Figure 3. Euclidean normalized difTerential number count comparison. 
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